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metrics 



o 
o 

(N 



o 
m 
o 

o 

7' 

u 

bJO. 



John D. Steele 

School of Mathematics, University of New South Wales, Sydney, New South Wales 2052, 
Australia, email: j.steele@unsw.edu.au 

Abstract: In this paper I show that the Newman- Tamburino spherical metrics always 
admit a Killing vector, correcting a claim by Collinson and French, (1967 J. Math. Phys. 
8 701) and also admit a homothety. A similar calculation is given for the limit of the 
Newman- Tamburino cylindrical metric. 

1 Introduction 

The Newman- Tamburino metrics are those vacuum solutions of the Einstein equations 
admitting hypersurface orthogonal geodesic rays with non- vanishing shear and divergence. 
In the Newman-Penrose formalism this implies that \E'o = k = 0, that p is real and non- 
zero and (T 7^ 0. In [1] Newman and Tamburino explicitly gave all such metrics and 
showed that they fall into two classes: the spherical, with ^ aa and the cylindical with 
= aa. In [2] Collinson and French claimed to have shown that the former metrics admit 
at most one Killing vector, and that happens only in a particular subcase. In fact, the 
spherical Newman- Tamburino metrics always admit a Killing vector and also always admit 
a homothety. This preprint is intended to show the full calculations and results when the 
homothetic equations of [3] are integrated for the Newman- Tamburino spherical metrics. 
The bulk of sections |21 and 0] come from Maple 9 worksheets, exported to T[t;X and suitably 
tidied up for better readabihty. 

Throughout I use the spin coefficient notation of [4]. For example I use k', p', a' and 
r' in place of the more traditional — z/, —p, —it and —A. 



2 Results 



The contravariant form of the Newman-Tamburino spherical metric [1] (see also [5] , equa- 
tion (26.21)) is 
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Here our coordinates are = m, = r, + ix^ = x + iy = ( and 

A{u) = bu + c, L = -\og('-^) a = A{CCy/^ R' = r^-a^. 

2 \r — a J 

Here h and c are real constants. 

The Collinson and French result (also quoted in [5]) is that there is a Killing vector 
only in the case where A is constant — in this situation the Killing vector is the obvious 
du- However, if 6 7^ we can set c = by a coordinate change and then the vector 

= -udu + rdr + 2x8^ + 2ydy. 

is a Killing vector, as will be shown in section 121 This can be checked directly: consider 
the flow of K"^. This scales the coordinates by 

u X^^u, r ^ Xr, ^ — > A^C 

for real parameter A > 0. Under this scaling it is easy to check that all the contravariant 

components given above are homogeneous in A (when A = hu) , and all of the correct degree 

to make the flow isometric. For example, the g^"^ component is homogeneous of degree 2, 

d d 

and so the metric term g'^^^r- ® tt is unchanged under the flow. 

or or 

Also the vector 

H = rdr + + ydy 
is a homothety, whatever A is (see section E}. Alternatively, the flow of H is 

u u, r ^ Xr, ( ^ AC, 

and we again flnd that all the contravariant components given above are homogeneous in A, 

and all of the correct degree to make the flow homothetic. For example, the g"^"^ component 

d d 

is homogeneous of degree 1, and so the metric term g"^"^— (8> 7:— scales by A ^: the same 

or or 

scaling applies to all the metric terms. 

Newman and Tamburino [1] also give the following metric, which arises as a limit of 
the cylindrical case (see also [5] (26.23) for corrections to the dv? coefficient): 

ds^ = 2dudr — x'"^ \b + log(r^x^)] du^ + A—du dx — r'^dx^ — x^dy"^ , 

with the same coordinates as used in the spherical case. The Killing vectors here are 
obvious {du and dy) and as we shall see there is also a homothetic vector (see section 

H2 = 2rdr — xdx + 2ydy. 

One can use the flow of H2 to check it is a homothety as well. 



2 



3 The Calculations (spherical case) 



The basic information is taken from Collinson and French [2], and Newman and Tamburino 
[1]. See those papers for those spin coefficients that are not actuaUy calculated here. I 
have checked in a separate calculation that their results are correct as quoted. I use as 
coordinates u, r, ( = x + iy. 

Collinson and French [2] wrote the conformal Killing equations in Newman-Penrose 
form and used that in their work, although there are a few minor typos in their paper. 
Here, I will use the formalism of [3], which generalised the ideas of [6] into a form suitable 
for this task. I will use to the notation of [3] for the components of the homothety 

and its bivector. Fab, with anti-self dual 

"Fab = 2000 ^[a^6] + 2001 (^[a^ft] - "^[a^b]) - 20ii n[„m6]. 

The tetrad is a standard tetrad (see [1]), based around the Debever- Penrose vector — dr, 
see [1] and [2] for further detail. Since the tetrad is normalised, for the Penrose-Rindler 
spin coefficients used in [3] we have 7' = — e, (3' = —a etc. 

In the Maple I use use z for ( and w for C; HI for ^£ etc. I typically add a b for a 
complex conjugate (Hmb is ^rn) ^-nd a 1 for a dash (rhol is p'). 

Firstly, define the terms a, aO (the latter is cco in [2]). 

> a:=A(u)*z^(l/2)*w~(l/2) : 

> aO : =3/4*w^ (3/4) (-1/4) : 

> aOb : =3/4*z^ (3/4) *w" (-1/4) : 

Rather than use the explicit definitions for L and R in [2] and [1] , 1 will leave them as 
"unknown" functions and define a routine later that will substitute for their derivatives. 
I will also use Q{u,r, z,w) in place of to make things more transparent. I define 
what these functions actually are so we can substitute for them more easily when that 
become useful. I also define dummy symbols to use in place of the full functional 
dependence of L and Q for ease of readability. 1 have also suppressed the functional 
dependence in the Maple output, replacing Q{u,r, z.iv) with Qlx"") for example. 

> LL:=L(u,r,z,w) :Lis:=l/2*log((r+a)/(r-a)) ; 

2 \r-A{u)^^J 

> QQ:=Q(u,r,z,w) : Qis : =1/ (r~2-a"2) ; 

Qis := (r'^ — {A (u))'^ zw) 
Now we define the routine to simplify derivatives and products and also add a line 
to collect terms. 
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> diffsbs:=proc(XX) 

> subsCdif f (L(u,r,z,w) ,r)=-a*Q(u,r,z,w) ,XX) : 

> subsCdiff (L(u,r,z,w) ,u)=r*Q(u,r,z,w)*dif f (a,u) ,%) ; 

> subs(diff (L(u,r,z,w) ,w)=r*Q(u,r,z,w)*diff (a,w) ,%) : 

> subsCdiff (L(u,r,z,w) ,z)=r*Q(u,r,z,w)*diff (a,z) //„) : 

> subsCdiff (Q(u,r,z,w) ,r)=-2*r*QCu,r,z,w)"2,yo) : 

> subsCdiff CQCu,r,z,w) ,z)=2*QCu,r ,z,w) "2*a*dif f Ca,z) //„) : 

> subsCdiff CQCu,r ,z,w) ,u)=2*QCu,r ,z,w) "2*a*dif f Ca,u) //o) : 

> subsCdif f CQCu,r,z,w) ,w)=2*QCu,r,z,w)''2*a*diff Ca,w) //o) : 

> student [powsubs] Cr'"2=a"2+l/QCu,r ,z,w) ,expandC/o)) ; 

> collect C/o , [L Cu , r , z , w) , Q Cu , r , z , w) , r , HI Cu) , psi , z , w] ) ; 

> end proc: 

The terms S and Sb are ip^ and its conjugate in [2]. 

> S : =2*A Cu) "2*z'^ C3/4) *\}^ C3/4) *z ; Sb : =2*A Cu) "2*z'' C3/4) kw" C3/4) *w : 

S := 2A{ufz''/^w^/^ 
And the curvature component ^'i is given in [1]. 

> Psil:=S*QQ''2; 

Now from [2] we have K = e = r' = \l/o = 0, and p and a real — these can also be easily 
checked by Maple. So by [3] (6a) D^e = 0. Using r = a + /3 ([2]), [3] (6c) becomes 

SCe = T^e - pirn - Cr^m " 011- 

We next use equation [3] (11), since is a Debever-Penrosc vector. Unfortunately, [3] (11) 
contains an error — the right hand side is the complex conjugate of what it ought to be. 
With this correction, we have 

011 = -Tie + pirn + (yim- 

Hence 5E,t = and i^ = ie{u), as found in [2]. 
Equation (10a) of [3] is 

D0n = -6^1, 

and so integrates to give the r dependence of (pu, here called pll. We ignore the factor 
independent of r when integrating: 

> intC-Qis"2,r) ; 



+ ^arctanh ( . , ^ ) A {u) ^ z ^ 



The arctanh term here is just L and we get 

> pll:=S*HlCu)/2/a"2*r*QQ-S*HlCu)/2/a~3*LL+pllOCu,z,w) ; 

pll := ^ + pllO{u,z,w) 

w'-'^ A [u) w^l^ 
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Here pllO(u,z,w) is the integration constant. Now the spin coefficients — see [2]. 

> alpha: =expand(simplify(S*LL*r/2/a'^2*CjQ+QQ*(r*a0+a*a0b)-S*QQ/a/2, radical)) : 

> beta : =-S*LL*QQ/2/a-QQ* (r*aOb+a*aO) : 

Now a routine to take conjugates nicely, as we need conjugates to define r. 

> conj :=proc(XX) 

> subs(z=wl,w=zl,XX) : 

> subs(wl=w,zl=z,y„) : subs(I=-I,y„) ; 

> subs (L (u , r , w , z) =L (u , r , z , w) , Q (u , r , w , z) =Q (u , r , z , w) , %) ; 

> subs(pllO(u,w,z)=pllOb(u,z,w) ,pllOb(u,w,z)=pllO(u,z,w) ,%) ; 

> end proc: 

> conj (alpha) +beta: 

> tau:=collect (%, [L(u,r,z,w) ,Q(u,r,z,w)] ) ; 

T := {-A (u) ^5/4^^1/4 + ^^3/4^-1/4-) Q (^a) ^ (^a) _ ^ ^i;5/4^1/4g (^a) 

> rho : =-r*Q (u , r , z , w) : 

> sigma : =a*Q (u , r , z , w) : 

> sigmal:=r*LL^2*S^2/4/a"4*QQ+r*LL*S*Sb/2/a"4*QQ-r*diff (a,u)*QQ-S*Sb/2/a'"3*CjQ: 

> unprotect (gamma) ; 

> gamma: =-r*LL"2*S*Sb/4/a~4*QQ+S~2*LL~2/4/a~3*QQ+a*r*LL*(S*a0b-Sb*a0)/2/a"3*QQ 

> + S*Sb*QQ*LL/2/a~3 - S*Sb/4/a^3*(LL/2/a~2-r/2/a*QQ) + 

> (S"2/2/a-a* (S*aOb-Sb*aO) ) /2/a~2*Qq : 

We need derivative operators 5 and 5' to find p'. Firstly, the components of 
come from [2] and [1]. 

> omO:=-A(u)*z~(l/4)*w"(5/4) : 

> omega:=-Sb*LL/2/a"2+(r*omO-a*conj (omO))*QQ; 

uj ■= -nfil^z-^I^L (x'^) + [-rA {u) z^'^w^'^ + A {uf z^l^nfil") Q (x") 

> omegal :=conj (omega) : 

> P:=z"(3/4)*w"(3/4) : 

> del : =XX->omega*dif f (XX , r) +2*r*P*Qq*dif f (XX , w) -2*a*P*QQ*dif f (XX , z) : 

> dell : =XX->omegal*dif f (XX , r) +2*r*P*QQ*dif f (XX , z) -2*a*P*QQ*dif f (XX , w) : 
To calculate p' we use [4] (4.11.12 e'). 

> expand((diff(sigmal,r)-rho*sigmal) /sigma) : 

> rhol:=diffsbs(y.) ; 

pi := -A (u) z^L {x^^f Q {x") -2iA{u)zw+ Q (x") L(x") +A{u)zw — Q {x") 

\ ^Jw J du 

We check some curvature equations next before we go on. 

> diffsbs(diff (tau,r) - rho*tau-sigma*conj (tau)-Psil) ; #[4]4.11.12c 
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> diffsbsCdiff (alpha, r) -rho*alpha-sigina*beta) ; # [4] 4.11.12h&i'; 



> diffsbsCdiff (beta,r)-sigina*alpha-beta*rho-Psil) ; # [4] 4. 11 . 12h'& i; 



> Psi2:=-diff sbsCdiff (rhol,r)-rhol*rlio-sigma*sgmal) ; # [4] 4.11.12 f 

^2 := -4A {uf z^'^L Q {x'^f ^-{2 A {u) z'^r + AA {uf z^l^w^/^) Q (x")' 
This expression for \I^2 agrees with [1]. 

> diffsbsCdiff (gamma, r) -beta*conj (tau)-alpha*tau-Psi2) ; # [4]4.11.12k 



> dif f sbs (del (rho) -dell (sigma) -rho* (conj (alpha) +beta) + 

> Sigma* (3*alpha-conj (beta) ) + Psil) ; # [4] 4.11.12 d 



> dif f sbs ( dell (beta) -del (alpha) -rho*rhol+sigma*sigmal+alpha*corLj (alpha) + 

> beta*conj (beta)-2*alpha*beta-Psi2 ); # [4] 4.11.12 1 



Integrating [3] {6g) and using [3](11) (corrected, see above): 

> Hml:=-r*pllO(u,z,w) + Hl(u)*S/2/a~3*r*LL+HmO; 

, HI (u) z^/^rL (x^) 
Hml := -rpllO (u, z, w) + — ^ , ' + HmO 

> Hmbl:=-r*pllOb(u,z,w) + Hl(u)*Sb/2/a'^3*r*LL+HmbO; 

HI (u) w^^^rL (x"") 
Hmbl := -rpllOb (u, z,w) ^ ' + HmbO 

By [3] (11) the following ought to be zero. 

> dif f sbs(pll+tau*Hl(u)-rho*Hml-sigma*Hmbl) : 

> collect ("/o/QQ, r) ; 

/ z^/'^Hl (u ) \ 
{HmO-\ ^^^^ ^ +A{u) ^/z^/wpllOh{u,z,w)\ r 

- A {u) y^VwHmbO - {A {u) f zw pllO {u, z, w) - HI (u) A (u) w''^^z^/^ 

> expcind(solve(coef f (yo,r , 1) ,HmO) ) , expand (solve (coeff (yo,r ,0) ,HmbO) ) ; 

z^/'^Hl (u) w'^/^Hl iu) 

— A iu) yGy/wpllOb {u, Z, w) , —y/zy/wA (m) pllO {u, Z, w) 



So we get 

> Hm : =-r*pl 10 (u , z , w) -a*pl 10b (u , z , w) +H1 (u) *expand (S/2/a~3* (r*LL-a) ) 
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> Hmb : =conj (Hm) : 

These agree with the components in [2] (their Vs and V4). Now we use [3] (10b) to 
get 001- 

> dif f sbs(Psil*Hm/2/sigma-beta*pll/sigma+del(pll)/2/sigma) : 

> p01:=collect(yo, [Hl(u) ,Q(u,r,z,w) ,L(u,r,z,w) ,r] ) ; 

p« := ( (i ^VSr + (2 + A (.) .») L Q - - ^] HI (,.) 

[U, Z, W ) 



+ 



A {u) z^l^w^l\ pi 10 {u, z,w)+A {uf z^/^w^/^pllOb {u, z, w)^ Q (x") 

/ 3pll0{u,z,w) wy'i^pllO{u,z,w) \ I _ ^3/4^3/4 A . ^ . 
\ ^A{u)w^l^ ^ A{u) y r z w ^^pnu[u,z,w) 

Now [3] (10c) and (8a) will give us information on the w (that is, Q dependence of 

> diffsbs(Psil*Hmb-Psi2*Hl(u)+2*rho*p01+2*alpha*pll-dell(pll)) ; # [3] 10c 

f SpllO {u,z,w) w^/^ ^ iin( \\ 1/4 
n—^ — TTT-'- — 2 ———-—pllO [u, z.w) \ z ' 

> dsolve(yo=0,pllO(u,z,w)) ; 

iin( \ -FUu,z) 
pllu [u, Z, W) - 



> diffsbsCdiff (p01,r)+dell(pll)-2*rho*p01-2*alpha*pll) ; # [3] (8a) 

f 9pllO{u,z,w) , ^ w'/^ d \ 
, . , , — T72 1- 3 ———pllO {u, z, w) z^i^ 

> dsolve(%=0,pllO(u,z,w)) ; 

-Fi {u, z) 
pllO [u, z, w) = ^ — 

Both giving the same result. Now we turn to ^„ and [3](6i), which we solve for cr^^. 

> rhs6i : =dif f sbs ( (-del (Hm) -conj (sigmal) *H1 (u) -Hm* (conj (alpha) -beta) ) ) 
The imaginary part ought to be zero as is real, so using results from [3] (10c) and 

(8a), we find the imaginary part divide out a common nor-zero factor and call what's 
left X. 



7 



> Imrhs 6 i : =%- c on j ( % ) : 

> subs (pi 10 (u , z , w) =F (u , z) /w^ (3/4) , pi 10b (u , z , w) =Fb (u , w) /z~ (3/4) , %) 

> X : =expand C/o/r/QQ/sqrt (z) /sqrt (w) ) ; 



X := -iz^^^Aiu) —F (u,z)+4:W^^^A(u) —Fb (u,w)+3 , -3 , 

oz aw w^/^ z^l^ 

Assuming F is differentiable in z we can split this 

> subs(Fb=0,X) :%; 

-4 z^l^A {u) —F {u, z) - 3 

This is a (real) function of u and w. We choose the shape of the separation function 
to simplify the solution to the differential equation slightly. 

> dsolve (%=-4*A (u) *G (u) , F (u , z) ) ; 

Check this out: 

> subs(F(u,z)=G(u)*z"(l/4)+H(u)/z"(3/4) ,Fb(u,w)=G(u) *w" (l/4)+Hb(u)/w- (3/4) ,X) : 

> expand (%) ; 



So we define a simplification routine for 0°^. 

> P110sbsl:=proc(XX) 

> subs(pllO(u,z,w)=F(u,z)/w"(3/4) ,pll0b(u,z,w)=Fb(u,w)/z"(3/4) ,XX) ; 

> subs (F (u , z) =G (u) *z~ (1/4) +H (u) /-zT (3/4) , Fb (u , w) =G (u) (1/4) +Hb (u) N" (3/4) , %) ; 

> expand (%); 

> end proc: 

And check it works 

> P110sbsl(Imrhs6i) ; 



Turning to [3] (66) next, 

> eqn6b : =H1 (u) * (gamma+conj (gamma) )-conj (tau)*Hm-tau*conj (Hm)-pOl-conj (p01)+psi: 

> PllOsbsKdiffsbsC/.)) ; 

^ ^ > 2z 2w 

This ought to be a function of u only, so if = and we define a new simplification 

routine and test it out: 

> P110sbs2:=proc(XX) ; 

> expand( subs( pll0(u,z,w)= (psi-diff (Hl(u) ,u))*z~(l/4)/w~(3/4) , 

> pllOb (u , z , w) = (psi-dif f (HI (u) , u) ) (1/4) /z~ (3/4) , XX) ) ; 

> collect(°/o, [L(u,r,z,w) ,Q(u,r,z,w) ,r,Hl(u) ,z,w] ) ;end proc: 

> P110sbs2(diffsbs(eqn6b)) ; 
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This is as it should be. Now we can define 

> P110sbs2(diffsbs(rhs6i+conj (rhs6i) )/sigma/2) : 

> Hn:=collect(yo, [L(u,r,z,w) ,Q(u,r,z,w) ,Hl(u) ,diff (Hl(u) ,u) ,r,psi] ) : 

We check this against the [2] version, called V2 there. It is clear from the shape of 
Cm {=V3 of [2]) that flo in [2] is my 0°^. 

> ayO:=pllO(u,z,w) :ayOb:=pllOb(u,z,w) : 

> V2:=r*LL"2*(-S"2-Sb"2)*Hl(u)/4/a~5 - LL"2*S*Sb*Hl(u)/4/a"4 - 

> LL*(ay0b*S+ay0*Sb)/2/a + r*(2*a*Hl(u)*diff (a,u)-ay0*Sb-ay0b*S)/2/a~2 + 

> (-2*a~3*ay0*S-2*a~3*ay0b*Sb+Hl(u)*S*Sb)/4/a"4+r*LL*(-2*a~3*ay0*S- 

> 2*a~3*ay0b*Sb- Hl(u)*S*Sb)/4/a~5 + l/QQ*(-3*ay0*S-3*ay0b*Sb- 

> 4*a^2*2*P*(diff (pllO(u,z,w),w) +dif f (pllOb(u,z,w) ,z) ) )/8/a~3: 

> expand(Hn-P110sbs2(V2)) : 

> simplify(subs(psi=0,diff (Hl(u) ,u)=0//.)) ; 



So our Cn agrees with [2] in the case of their Killing vector (-0 = and $,e constant) . 
However, if C,e is not constant, the terms differ: 

> simplify(subs(psi=0,°/o°/o)) ; 



Next, we put our ^„ into [3](6(i). 

> eqn6d:=diff sbsCdiff (Hn,r)-p01-conj (pOl)-psi) : 

> P110sbs2(%); 

A(u) du ^ ' 



> dsolve(y.,Hl(u)); 

HI (u) = _C1 A {u) 

So next a routine to replace Ci{u) with a multiple of A{u), and also to kill off the 
second derivative of A[u). 

> Hlsbs:=proc(XX) 

> subs (HI (u) =C*A (u) , XX) ; subs (dif f (A (u) , u , u) =0 , %) ; subs (dif f (HI (u) , u , u) =0 , %) ; % ; 

> end proc: 

We now try [3](6j). 

> dell(Hm)+conj (rliol)*Hl(u)+rho*Hn+(conj (beta) -alpha) *Hm-p01+conj (p01)+psi: 

> Hlsbs (PI 10sbs2(diffsbs (%))); 
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So that is satisfied. Now for 0oo, which we get from the conjugate of [3] (6/). 

> eqn6f :=-dell(Hn)-(conj (beta)+alpha)*Hn-conj (rhol)*Hmb-sigmal*Hm: 

> p00:=diffsbs(P110sbs2(diffsbs(%))) : 

I've supressed this component as it's very long, but we check the result with [3] 

(Sd). 

> dell(p01)+diff (pOO,r)-rho*pOO-sigmal*pll: 

> P110sbs2(diffsbs(yo)) ; 



To go any further we need to get the components of n" = {1,U,X^,X^) and to 
define D'. Taking the metric terms from [1] and [2]: 

> gup22:=-2*r~2*sqrt(w)*sqrt(z)*QQ + 2*r*LL/A(u) + QQ'^2* (2*r'^3*A(u) * (w"2+z~2) 
- 4*r"2*A(u)"2*w"(3/2)*z"(3/2)) : 

> gup22+2*oinega*conj (omega) : 

> U:=diffsbs(y./2); 

U : = V^V^L {x'^f +{{w'' + z"") A{u)r-2A {uf z^/^w^/^) Q - V^V^ 
+ (^{-A [uf z-'/^v^ + 2 zwrA (u) - A {uf ^zvS'l'^) Q {x") + JJ^^ L (x") 

> gupSa : =-2*z~ (3/2) *w" (3/2) / (r+a) "2 : gup44 : =-2*z" (3/2) *w~ (3/2) / (r-a) "2 : 

These next two terms are the components of m". 

> xi3 : =P* (r-a) *QQ ; xi4 : =I*P* (r+a) *QQ ; 

Xi3 := ^3/4^3/4 _ ^ (^^^ y/Zy/w) Q (x") 

xi4 := iz^/^w^/^ {r + A (u) y/z^) Q (x") 

> xi3*conj (xi4)+xi4*conj (xi3) ; # checking 



> simplify(( subs( Q(u,r,z,w)=Qis,xi3*conj (xi3)*2+gup33) )); 



> simplify (subs (Q(u,r,z,w)=Qis,xi4*conj (xi4)*2+gup44)) ; 



> gup23 : =4*A(u) ~2*z- (3/2) *w" (3/2) * (z+w) /2* (LL/2/a"3- (r-2*a) *QQ/2/a"2 

> -(r-a)*QQ~2) : 

> gup24 : =4*A(u) "2*z^ (3/2) (3/2) * (z-w) /2/I* (LL/2/a''3- (r+2*a) *QQ/2/a"2 

> -(r+a)*qQ~2) : 

> gup23+omega*conj (xi3) +conj (omega) *xi3 : 

> X3:=diffsbs(y.) ; 
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> gup24+omega*conj (xi4)+omegal*xi4: 

> X4 : =dif f sbs (factor (%) ) ; 

X4 := i |^((-^3/2^ ^ r + zw^A (u) - z^wA (u)) Q + - L (x'^) 

+ + iz^^'^y/w) r + i^^wA (u) - izw'^A (u)) Q {x") 

As a double check we firstly define the (contravariant) tetrad and then check against 
the metric terms. 

> ell : =<0 , 1 , , 0> : en : =<1 , U , X3 , X4> : 

> em :=<0, omega, xi3,xi4>:emb:=map( con j ,em) : 

> ell . Transpose (en) -em . Transpose (emb) : 

> %+Transpose(yo) : 

> g:=map(diff sbs,%) : 

> diffsbs(g[2,2]-gup22) ; 



> diffsbs(g[2,3]-gup23) ; 

> diffsbs(g[2,4]-gup24) ; 










> diffsbs(simplify(g[3,3]-gup33)) ; 



> diffsbs(simplify(g[4,4]-gup44)) ; 



For a second check we apply two of the commutators [4] (4.11.11) to r and check 
what we get. 

> dif f (U,r)+gamma+conj (gamma) -tau*conj (omega) -conj (tau)*omega: 

> dif f sbs (7.) ; 



> dif f (X3 , r) -t au*conj (xi3) -conj (tau) *xi3 : 

> dif f sbs (7.) ; 
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Since all this checks out we go ahead and define D' . 

> Dl:=proc(XX) 

> dif f (XX , u) +dif f (XX , r) *U+ (X3+I*X4) *dif f (XX , z) + (X3-I*X4) *dif f (XX , w) ; 

> P110sbs2(diffsbs(%)) ; 

> end proc: 

Wc make use of D' firstly to find the last spin coefficient, k', using [4](4.11.125f). 

> Dl (beta) -del (gamma) -tau*rhol-alpha*conj (sigmal) -beta* (rhol+gamma-conj (gamma) ) 

> +gamma*(tau-beta-conj (alpha)) : # should be -kappa' *sigma 

> diffsbsC/.): 

> kappal:=-diffsbs(P110sbs2(y./sigma)) ; 



-2^ + 2 z'/^w'/^^ r - 3 z^^w^'^A {u)^ Q L + 
1^^1/4^5/4 _ 2 ^5/4^1/4^ r-A{u) (^^7/4^3/4^ ^ 4z'/^w'/^^ ^ Q (x") L (x") 

- (-2 z'/^w^/'A (u) + w'/^^/'A {u)^-2 zV4^l/4^1jl^ g ^^a^ 
\ du du J 

Now to look at the equations that involve D' . Firstly [3] (6/i): 

> eqn6h:=Dl(Hm)+conj (kappal)*Hl(u)+tau*Hn+(conj (gamma) -gamma) *Hm- con j (pOO) 

> P110sbs2(diffsbs(%)) :factor(Hlsbs(y.)) ; 



Then we look at [3] (6e): 

> eqn6e : =D1 (Hn) + (gamma+conj (gamma) ) *Hn+kappal*Hm+conj (kappal) *Hmb : 

> P110sbs2(diffsbs(%)) : 

> factor (Hlsbs (%)) ; 



and [3] (8c) 

> eqSc : =dif f sbs (del (pOl) +D1 (pll) -sigma*p00-2*tau*p01-(rhol+2*gamma) *pll) : 

> factor (Hlsbs (P110sbs2 (%) ) ) ; 



and [3] (86) 

> Dl (p01)+del (pOO) - (tau-2*beta) *p00-2*rhol*p01-kappal*pll : 

> P110sbs2(diffsbs(y.)) : 

> factor (Hlsbs (%)) ; 



and finally [3] (lOd). 
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> eqnl0d:=Psi2*Hm-Psil*Hn-2*tau*p01-2*gainma*pll+Dl(pll) : 

> P110sbs2(diffsbs(%)) :factor(Hlsbs(y.)) ; 



Next we consider what happens if we have a KiUing vector (-0 = 0) with zero 

iC = 0) 

> subs(C=0,psi=0,Hlsbs(P110sbs2(Hm))) ; 



> subs(C=0,psi=0,Hlsbs(P110sbs2(Hn))) ; 



Hence we cannot have both ip and C zero. This is the CoUinson and French result: 
only one Killing vector at most. We have a look at the homothety. 

> HI (u) *en+Hn*ell-Hin*emb-conj (Hm) *ein : 

> mapCdif f sbs//o) : 

> map(P110sbs2,yo) : 

> map(Hlsbs,%) : 

> subs(z=x+I*y,w=x-I*y,K) :K:=map( expand//,) ; 



K :-- 



( CA{u) \ 

-rC^ + 2r^ 

^ an 



So the obvious Killing vector if A is constant: 

> KK:=subs(psi=0,C=l/B,diff (A(u) ,u)=0,A(u)=B,K) :%; 

/1\ 





voy 

The new Killing vector in the other case: 

> KK2:=subs(psi=0,A(u)=u*B,C=l/B,K) : map (simplify, 7.) ; 

( u\ 

—r 
-2x 

And the proper homothety for both cases: 

> HH:=subs(C=0,psi=l,K) ; 
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2r 
2x 



We now turn to the remaining curvature equations and Bianchi identities. To make 
life easy, we define the weighted derivative operators, [4] section 4.14. 

> thorn:=X->diffsbs(diff (X,r)) : 

> thornl :=proc(X,p,q) 

> local i; 

> Dl(X)-p*gaiifflia*X-q*conj (gamma) *X; 

> dif f sbs (expand(yo) ) ; 

> end proc: 

> edth:=proc(X,p,q) 

> local i; 

> del (X) -p*beta*X-q*conj (alpha) *X ; 

> diff sbs (expand (%)) 

> end proc: 

> edthl:=proc(X,p,q) 

> local i; 

> dell(X)-p*alpha*X-q*conj (beta)*X; 

> diff sbs (expand (%)) 

> end proc: 

As a check on the calculations, we can run through the curvature equations, [4] 
(4.12.32), some of which we've used already, some of which will give us ^3 and ^4. The 
only ones that do not give zero are (6') and (c'), the first of which gives us ^4: 

> thornl (sigmal , -3, 1) -edthl (kappal ,-3,-l)-sigmal*(rhol+conj (rhol) ) 

> +kappal* (conj (tau) ) : 

> Psi4:=Hlsbs(diffsbs(y.)) ; 




du 



And (c') gives '^^ (as do several others): 
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> -thorn (kappal)+conj (tau)*rhol+siginal*tau: 

> Psi3:=diffsbs(y.) ; 



+ (^14 {A {u)f z^l^uf^l^ + 6 ^M^^^ g [x'^f + 2 Q z'^l^w'l^ L 
+ (2 A {uf z'^/^w'/^^ + 6 z'/^w'/^A {uf + 6A (u) z'/^w^^^r^ Q {x'^f 

+ 2^z'/W/^Q{x^) 
du 

Now we check the leading terms (in inverse powers of r) of our and and 
compare to [1] . 

> subs(L(u,r,z,w)=Lis,Q(u,r,z,w)=Qis,Hlsbs(Psi4)) : 

> T4:=subs(r=l/R,yo) : 

> series (T4,R=0, 3) assuming R: : positive; 

-8^zwR^ + 0(R') 
du ^ ' 

Here the leading term agrees with [1] . Next 

> subs(L(u,r,z,w)=Lis,Q(u,r,z,w)=Qis,Hlsbs(Psi3)) : 

> T3:=subs(r=l/R,%) : 

> series (T3,R=0, 4) assuming R: : positive; 

2 1^^5/4^1/4^2 ^ 8^7/4^3/4^ ^3 ^ ^ /^4) 

du ^ ' 

We find that the leading term agrees with [1], but in the second term the powers of z = C 
and w = are wrong in [1]. We can also check that the Bianchi identities, [4] (4.12.36-39) 
are satisfied (and they are). 

Finally, we turn to the remaining integrability conditions, [3] (lOe) to (lO/i). 

> P110sbs2(diffsbs(Psi2*Hmb-diff (p00,r)-Psi3*Hl(u))) ; # [3] lOe 



> Psi3*Hm-Psi2*Hn-2*rhol*p01+2*beta*p00+del(p00) : # [3] lOf 

> P110sbs2(diffsbs(%)) ; 



> Psi3*Hmb-Psi4*Hl(u)+2*sigmal*p01-2*alpha*p00-dell(p00) : # [3] lOg 

> P110sbs2(diffsbs(%)) ; 



> Psi4*Hm-Psi3*Hn-2*kappal*p01+2*gamma*p00+Dl(p00) : ## [3] lOh 
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> Hlsbs(P110sbs2(diffsbs(y.))) ; 



So we see that all the homothetic and Killing equations are satisified and we have shown 
that there is always a Killing vector in these metrics and also always a homothety. 

4 The Calculations (limit cylindrical case) 

Since neither [1] not [2] give the spin coefficients for the limit cylindrical metric, we will 
need to calculate them using Maple's tensor package. Note that we use the corrected 
version of this metric, see [5] equation (26.23) 

> with(tensor) : 

> coord: = [u,r ,x,y] :g_c : =array (1 . .4,1. . 4, symmetric , sparse) : 

> g_c[l,l] :=-expand(simplify((b+log(r~2*x"4))/x"2/2) 

> assuming r : :positive,x: :positive) ; 



> g_c[l,2] :=l:g_c[3,3] : =-2*r-2 : g_c [4,4] : =-2*x"2 : g_c [1 , 3] :=2*r/x: 

> g:=create( [-1,-1] ,eval(g_c)) : 
Next we calulate all the relevent tensors. 

> tensorsGR(coord,g,gup, 'detg' , 'CI' , 'C2' , 'Rm' , 'Rc' , 'R' , 'G' , 'C ) ; 
To calculate the spin coefficients, I use a set of routines available on my web site 

http : / / www . maths . unsw . edu. au/f ^j ds/papers . html 

> read 'PRcoeff ' : 

Now we define the tetrad, with the choice of dictated by the need for the tetrad 
to be right-handed, so the anti-self duality used in the definition of the homothetic 
bivector (see [3]) is satisfied. 

> md:=create( [-1] ,vector( [0,0,r,-I*x] )) :mup:=raise(gup,md, 1) : 

> mbd: =create( [-1] , vector ( [0 ,0,r , I*x] ) ) :mbup : =raise(gup,mbd, 1) : 

> ld:=create( [-1] , vector ( [1,0,0,0])) : lup : =raise (gup, Id, 1) : 

> nd:=create( [-1] , vector ( [g_c [1 , 1] /2, 1 , 2*r/x,0] ) ) :nup: =raise(gup,nd, 1) : 

> his :=linalg[stackmatrix] (ld[compts] ,nd[compts] ,md[compts] ,mbd[compts] ) : 

> h:=create( [1,-1] ,op(his)) : 

Using the routines PRspin and PRcrv from the PRcoef f file we calulate the spin 
coefficients, and curvature components. 

> spins :=PRspin(g,h,C2, coord) : 

> crv :=PRcurve(g,h,C,Rc, coord) ; 



b 
2^ 



In (r) 
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Prom these two calculations we find that the non-zero spin coefficients are 



1 1 1 , , 6 + log(rV) 

/3 = T=- — , p = (7=- — , 7=-- ^, p=(7 = 



2rx' 2r' 4ra;^' 8rx^ ' 

and the non-zero curvature components are 

1 ^ 1 ^ h + log(r^x^) 

^1 = ; ^2 = ; ^3 = ^ ■ 

Now we define the derivative operators D, D' , 5 and 5' 

> De:=XX->add(lup[compts] [i] *diff (XX, coord [i] ), i=l . .4): 

> Dl : =XX->add (nup [compt s] [i] *dif f (XX , coord [i] ) , i=l . . 4) : 

> del :=XX->add(mup [compt s] [i] *dif f (XX, coord[i] ) , i=l . .4) : 

> dell:=XX->add(mbup[compts] [i] *diff (XX, coord [i] ) ,i=l. .4) : 

Now to find the Killing vectors. Using [3] (6a), (6c) and (11) gives = ^e{u). Then 
from [3] (10a) we get 0ii, and find that 4>ii{u,x,y), the integration constant, is real by 
[3] (65^), which also gives So 

> Hin:=-r*pllO(u,x,y)+I*HmO(u,x,y) :Hmb:=-r*pllO(u,x,y)-I*HmO(u,x,y) : 

> pll:=Hl(u)/2/x/r+pllO(u,x,y) ; # note that pllO is real 

HI (u) 

pll ■.^^^ + pllO{u,x,y) 
We also solve [3] (106) for 0oi- 

> crv[Psil] *Hm-2*spins [sigma] *p01- (spins [beta] -spins [alphal] )*pll+del (pll) 

> pOl :=expand(solve(yo,p01) ) ; 



pllO{u,x,y) .HmO{u,x,y) Hl{u) Id .rd 

POI := ^ — h-—pllO{u,x,y)-i——pllO {u,x,y) 

2x 2xr Ax^r 2 ax 2x ay 

Now looking at [3] (8a), using the fact that k = 0: 

> dif f (pOl , r) +dell (pi 1) -2*spins [rho] *pO 1- 

> (spins [taul] +spins [alpha] -spins [betal] ) *pll : 

> expand (%) ; 

So 05i is independent of y. Now looking at [3] (66): 



> 



Dl (HI (u) )+2*spins [epsilonl] *H1 (u) +spins [tau] *Hm+spins [tau] *Hmb+p01 
+subs(I=-I ,p01)-psi: 

collect (°/o,r) ; 
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d d 

—HI iu) + ^V^IO {u, x,y)-iP 

So we solve this for recalling that 05i is independent of y, and use it to redefine 
)ii , 001 and ^rn- 

> pll:=Hl(u)/2/x/r+(psi-diff (Hl(u) ,u))*x+pO(u) : # note that pO is real 

> p01:=expand(subs(pllO(u,x,y)=(psi-diff (Hl(u) ,u))*x+pO(u) ,p01)) : 
Hin:=expand(subs(pllO(u,x,y)=(psi-dif f (Hl(u) ,u))*x+pO(u) ,Hm)) ; 



> 



Hm :— —rxtp + rx4-Hl [u) — rpO {u) + iHmO {u, x, y) 



' du 

> Hmb:=expand(subs(pllO(u,x,y)=(psi-dif f (Hl(u) ,u))*x+pO(u) ,Hmb)) : 

Turning to [3] (10c) 

> crv[Psi2] *Hl(u)-crv[Psil] *Hinb-2*spins [rho] *p01- (spins [alpha] -spins [betal] )*pll 

> +dell(pll): 

> expand (%) ; 



Now the right hand side of [3] (6d) is 

> -2*spins [epsilon] *Hn-spins [taul] *Hm-spins [taul] *Hmb+p01+subs(I=-I ,p01)+psi : 

> expand (%) ; 

du X 2x^r 

This is D^n, so we integrate 

> intC/o.r); 

^ UK \ n( \ Hl{u)\n {r) 

r—Hl [u) — 2 —pO [u) 



(III X 2x 



> Hn:=%+HnO(u,x,y) : 

Turning to [3] (6i) , 

> del(Hin)+spins [sigmal] *Hl(u)+spins [sigma] *Hn+ (spins [alphal] +spins [alpha] )*Hin: 

The coefficients of r are independent, so we collect the terms. 

> collect (expandC/o) ,r) ; 



-lb H -HI iu) - ^ ' —HmO iu, x, u) 

^ 2du ^ ' 2x 2xdy ^ ' '^^ 

fl 9 , Hliu)h Hliu)\Yiix) iHmOiu,x,y) 1„ , A , 

The imaginary part of the r"^ term implies = ^fiu, y), for some function f{u, y), 

so: 
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> X:=collect (subs(HmO(u,x,y)=x*f (u,y) ,%) ,r) 

> Y : =solve (op (5 , X) , HnO (u , x , y) ) ; 



Y := -HI {u) 



6 + 4 In {x) 



> expand (subs (HnO (u, X, y ) =Y, X) ) ; 



Id pO{u) Id 



> XX:=rhs(dsolve(%,f (u,y)))*x; 

XX := (-2 y^P + y^Hl {u) - + _F1 {u) \ x 

\ du X 

Where _F1 is an arbitrary function. This XX is So 

> Hin:=collect(subs(HmO(u,x,y)=XX,Hin) , [r,x,y] ) ; 

> Hinb:=collect (subs(HmO(u,x,y)=XX,Hinb) , [r,x,y] ) : 



Hm 



-ijj + —HI (u)]x- pO (u) 
du 



r+ 



M -2 ^ + —HI (u)]y + LFl (u) 
du 



x—iy pO (m) 



> Hn : =subs (HnO (u , x , y) =Y , Hn) ; 



du J X Ax^ 



2 In (r) + 6 + 4 In {x) 



> p01:=expand(subs(HmO(u,x,y)=XX,p01)) ; 

pO (u) iyib iy d , . iy , , i „ , , Hl(u) 
pOl := ^--^ + -^—Hl (u) + -^pO (u) - —_F1 iu) - —\^ 

2x r 2rdu ^ ' 2rx^ ^ ' 2r ^ ' Ax^r 

Returning to the integrabihty conditions, we look at [3] (lOd) 

> eqlOd:=crv [Psi2] *Hm-crv [Psil] *Hn-2*spins [tau] *p01-2*spins [gamma] *pll+Dl (pll) 

> expand (%) ; 

pO in) d ^ , , ^„ , , 

So by comparing coefficients we have 

> pO (u) : =0 ; HI : =x->kO*x+kl ; 

pO {u) := 

Hl:^ x^ kOx + kl 
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And a quick check shows that eqnlOd ([3] (lOd)) is satisfied. Next, the conjugate 
of [3] {6h) will give us 0oo- 

> Dl (Hmb) +spins [tau] *Hn+ (spins [gamma] +spins [epsilonl] ) *Hmb : 

> pOO:=collect (expand (%) , [psi,kl,kO]) ; 

00 ( ^ In (r) In ^ ^ /In (r) ^ In (x) ^ b \ 

\ 4:X 2x X J \ Arx^ 2rx^ 8rx^ J 

ub b 1 In (r) ln(x) u\n(x) ln(r)-u\ . f d , .\ 

H — H + . , ] kO - li —-F1 (u) X 



8rx^ Ax 2x 2x x 2rx^ Arx^ J \du 

We next check some further integr ability conditions, [3] (lOe) first. 

> crv [Psi3] *H1 (u) -crv [Psi2] *Hinb-2*spins [taul] *p01+2*spins [epsilon] *pOO 

> +diff (pOO,r) : 

> expand (%) ; 



And then [3] {8d). 

> e8d:=dif f (pOO,r)+dell (pOl) -spins [rho] *p00-2*spins [taul] *p01-spins [sigmal] *pll ; 

> expand (%) ; 

i—-Fl (u 

2r \du ^ ' 

So the intcgability function _F\ is constant: 

> _Fl(u) :=k3; expand (e8d) ; 

_Fi {u) := k3 


Also [3] (6e) is 

> eqn6e :=expand(Dl(Hn)+2*spins [gamma] *Hn) ; 

kO 

eqnoe 



2x2 

Thus kO = 0, and the coefficients simplify as follows: 
> Hl(u);Hn;Hm; 



kl 

ln(r) kl kl (6 + 4 ln(a;)) 
2^2 4a;2 
—rxip + (—2 iyip + ikS) x 
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> pOO;p01;pll; 



Ax 2x 



b In (r) 



In (x) 



X 




In (r) In [x) 



rx^ 2rx^ 



+ 



8rx^ 



b 



) 



kl 



iyip ikS 



kl 



r 2r 



Ax'^r 




All the remaining homothetic equations and integrability equations are satisfied, 

and we are left with the general homothetic vector: 

> lin_com(Hl(u) ,nup,Hn,lup,-Hin,nibup,-Hmb,mup) ; 

TABLE {[index _char = [Vl^compts = vector {[kl ,2rijj, —xijj,2yijj — k3])]) 
That is, 
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kidu + k^dy + {2rdr - xd^ + 2ydy) . 
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